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In  1968  Sherbrooke  formulated  the  well  known  METRIC  model  for  a two- 
echelon  inventory  system  consisting  of  a set  of  bases  and  a supporting  depot. 

The  items  stocked  in  the  system  are  called  recoverable  items,  that  is,  they 
are  subject  to  repair  when  they  fail.  Sherbrooke  and  others  have  proposed  a 
variety  of  algorithms  for  determining  optimal  base  and  depot  stock  levels. 

A substantial  portion  of  the  computational  requirement  associated  with  each 
of  the  algorithms  is  related  to  the  search  for  the  optimal  depot  stock  level. 

The  purpose  of  this  paper  is  to  describe  an  easily  implementable  method  for 
estimating  the  optimal  depot  stock  level.  The  computational  experience  re- 
ported in  the  paper  indicates  the  proposed  method  provides  an  excellent  estimate 
of  the  optimal  depot  stock  level  particularly  for  high  demand  items.  Further- 
more, the  proposed  method  significantly  reduces  the  computational  requirements 
for  any  known  algorithm  for  solving  Sherbrooke's  problem. 
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supply  organization  at  which  the. original  demand  was  placed.  All  resupply  is 
done  on  a one-for-one  basis. 

Sherbrooke  describes  both  the  problem  and  its  formulation  and  analysis  in 
detail  in  reference  2.  In  addition,  Sherbrooke  proposed  a marginal  analysis 
algorithm  for  solving  the  problem.  Another  approach  has  been  suggested  by  Fox 
and  Landi^  for  determining  optimal  base  and  depot  stock  levels.  A substantial 
portion  of  the  computational  requirement  associated  with  both  of  these  algorithms 
is  related  to  the  search  for  the  optimal  depot  stock  level.  As  will  be  discussed, 
this  search  is  particularly  time  consuming  for  items  having  a high  expected 
number  of  items  in  transit  to  the  depot  and  in  depot  repair. 

The  purpose  of  this  paper  is  to  develop  a method  for  estimating  the  optimal 
depot  stock  level  which  is  both  accurate  and  easy  to  implement.  Computational 
experience  reported  in  this  paper  indicates  the  proposed  method  provides  a very 
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good  estimate  of  the  optimal  depot  stock  level.  Furthermore,  the  proposed 

method  significantly  reduces  the  overall  computational  requirement  for  the 

• 1 2 
algorithms  developed  by  Fox  and  Landi  and  Sherbrooke  . 


2.  BACKGROUND 

The  optimization  problem  posed  in  reference  2 by  Sherbrooke  for  finding 
optimal  stock  levels  for  the  depot-base  system  described  in  Section  1 is 


'in  [ j I (y  ~ s..,)p(x|Ai1Tij(si0)) 


j=l  i=l  x>s. . 

m n 

subject  to  l c.  I s..  < C, 
i=l  j=0  13 

s • • — Ojij ii«) 

i: 


r'  3-v 
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where  n is  the  number  of  bases, 

m is  the  number  of  items, 

5^^  is  the  stock  level  at  location  j for  item  i, 

s^0  is  the  depot  stock  level  for  item  i, 

is  the  expected  daily  demand  rate  for  item  i at  base  j, 
c.  is  the  unit  cost  for  item  i, 

i 

C is  the  budget  constraint, 

Tij(Si0)  is  the  average  resupply  time  for  base  j for  item  i given 
the  depot  stock  level  for  item  i is  s^Q, 

-A..T.  (s.  ) 

and  p(x|X..T..(s.  ))  = e 1]  13  1U  (X..T..(s.  ))  'x!,  where  p(x|*) 

1]  lj  10  1]  iu 

represents  the  probability  of  having  x units  in  the  resupply 

system  at  any  point  in  time. 
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Furthermore,  Sherbrooke  shows  that  T. ,(s.„)  can  be  expressed  as 

W ■ v«  * a-rij,(!!ij  * s(sio)-Di>'  (2) 

where  A„  is  the  average  base  repair  time  for  item  i at  base  j, 

r„  is  the  proportion  of  demands  requiring  base  repair  for  item  i 

at  base  j , 

is  the  average  ordcr-and-ship  time  at  base  j for  item  i, 
is  the  average  depot  repair  cycle  time  for  item  i, 

A i _ 

6(si0)*Di  = ~ l (x-s^0)p(x|A^D^),  the  expected  depot  delay  per  demand 


for  item  i , 

A n 

A^  = l (l-r^)A^,  the  expected  daily  depot  demand  rate  for  item  i. 


Subsequently  Fox  and  Landi  proposed  a Lagrangian  formulation  of  Problem  1 
and  suggested  a relatively  simple  method  for  solving  the  Lagrangian  problem. 

This  approach  requires  only  a small  fraction  of  the  time  consumed  by  Sherbrooke's 
method  to  obtain  the  optimal  solution.  Their  Lagrangian  formulation  of  the 


problem  is 


u m u in 

III  (x-s  .)p(x]A  T (s.  ))  + e l l CS 
j=l  i=l  x>s . . 13  13  13  10  j=0  i=l  1 13 


subject  to  s^j  = 0,1,..., 

where  0 is  a Lagrangian  multiplier.  Fox  and  Landi  propose  to  solve  Problem  1 


I » •. 

!P 

t"*  ftp. 
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by  solving  Problem  3 for  a fixed  finite  set  of  multiplier  values  ©q,...,©^. 
In  their  approach,  the  ''optimal”  value  of  0 is  the  0 , Kc{0,...,M},  whose 

lx 

total  investment  cost  is  closest  to  C.  We  will  now  see  how  the  solution  to 
Problem  3 is  found  for  a given  value  of  0. 


4 


Observe  that  ProBlem  3 is  separable  by  item.  Consequently  for  each  6,  its 
solution  can  le  obtained  by  solving  m independent  item  problems.  To 
simplify  the  notation  we  drop  the  iter  designator  in  what  fallows.  Then  for 
each  item  for  a < given  9 we  solve  a problem  of  the  form 

n n 

min  l l (x-s . )p(x|X.T.(sn))  +6  l cs 


j=l  X>S_. 


1 rr  o 


subject  to  s ^ = 0,1,...,  where  j = 0, ,n.  This  problem  can  be  partitioned 

as  follows: 


n 

min  {0cs  + £ min{  £ (x-s . )p(x | X .T .(s  ) ) + 0cs . : s.s0,l,...,  s fixed}}. 

=0,1,...  j=l  s,  x>s.  3 3 3 33  (u) 


1 3 


Problem  4 is  not  necessarily  convex  in  sQ;  however,  the  inside  minimization 
problems  are  convex  for  a fixed  value  of  sQ.  Ke  may  restate  Problem  4 as 


follows : 


min  Z(so;0) 
Sg-0,1, . . . 


where 


Z(so;0)  = 0csQ  + ^ nun{J^  (x-sj}p(x|XjTj(s0))  + 0cSj:  Sj  = 0,1,...,  sQ  fixed}. 


3 3 


To  determine  Z(sQ;0),  we  must  solve  the  n base  problems  of  the  form 

min  ][  (x-s.)p(x|X.T.(s  ))  + 0cs..  (€ 

s . = 0,1, . . . x>s . 3 


It  is  easy  to  see  that  the  optimal  s^  for  Problem  6 is  the  smallest  nonnegative 
integer  for  which 


. IS? 

£ 


l p(x|X.T.(s  ))  < 0c. 

x>s.  3 3 
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To  find  the  optimal  solution  to  Problem  5 we  must  implicitly  examine 
Z(sQ;0)  for  sQ  = 0,1,...  . In  practice,  however,  the  number  of  values  for 
Sg  that  need  to  be  examined  explicitly  is  small.  This  is  due  to  the  asympto- 
tic behavior  of  6(sQ).  Specifically,  i$(Sq)  approaches  0 quite  rapidly 
as  sQ  increases  beyond  xD.  (Experience  gained  by  the  author  in  applying 
the  Fox-Landi  technique  to  real  Air  Force  problems  suggests  that  no  more  than 
15  values  for  sQ  need  to  be  explicitly  examined  for  any  item.)  Observe 
that  the  amount  of  time  required  to  solve  Problem  1 using  the  Fox-Landi 
algorithm  is  roughly  proportional  to  the  number  of  depot  stock  levels  for 
which  Z(sQ;0)  is  explicitly  evaluated.  Computation  time  required  by  Sher- 
brooke's algorithm  to  obtain  the  solution  to  Problem  1 is  also  proportional 
to  the  number  of  depot  stock  levels  examined.  Consequently,  if  the  number 
of  depot  stock  levels  that  are  examined  explicitly  can  be  reduced,  then  the 
total  time  required  to  find  the  optimal  solution  to  Problem  1 can  also  be 
reduced . 

In  the  next  section  we  develop  a method  that  provides  an  excellent  esti- 
mate of  the  optimal  depot  stock  level.  This  method  significantly  reduces  the 
computation  time  required  to  find  an  optimal  solution  to  Problem  1 using 
either  the  Fox-Landi  or  the  Sherbrooke  algorithm. 


3.  THE  ESTIMATION  PROCEDURE 

In  the  preceding  section  we  indicated  how  the  optimal  base  stock  level. 


call  it  sV,  can  be  determined  for  a given  value  of  sQ  (and  0).  In 
particular,  we  showed  that  sV  is  optimal  if  it  is  the  smallest  nonnegative 
integer  for  which 
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We  begin  this  section  by  developing  a different  but  equivalent  way  of  charac- 
tenzing  sj.  We  then  show  how  to  find  an  estimate  of  the  optimal  depot 
stock  level  by  solving  a relaxed  version  of  Problem  4 in  which  both  the 
nonnegativity  and  integrality  requirements  for  s^  are  dropped. 


Define  the  backorder  function  for  base  j for  a given  depot  stock  level 


so  as 


B.(s.;s  ) = ][  (x-s . )p(x| A .T .(s  )),  for  s.  t>  0 and  integer, 

j J V So  J ] 3 ^ 3 


and  the  convex  piecewise  linear  completion  of  B.,  call  it  B.,  as  follows: 

D j 

|^Bj(tjsQ),  if  t is  a nonnegative  integer. 


A J 

. (t;sn)  - / TB  . (s  . ;s  ) - B.(s,-l;i 


T ’ 0 


B^s^Sg)  - B j ( s -1 ; s0  )3  (t-  ( s^  -1)  ) + B^.(s].-l;s0),  s..-l<t<Sj; 
where  s^  is  a nonnegative  integer,  and  Bj(-l;sQ)  = °°. 


Let  ABj(Sj;sQ)  = B^s^.jSg)  - Bj(Sj-l;sQ)  when  s_.  is  a nonnegative  integer. 


D(Sj >sq)  = {v:  ABj(s^;s0)  < v 4 ABj(s_.+l;sQ)}. 


Observe  that  Dj(Sj;sQ)U{ABj(Sj',s0)}  is  the  set  of  subgradients  of  B^  at 
s..  Then  an  alternative  way  of  verifying  that  s.  is  an  optimal  base  stock 

A 

level  is  to  show  that  -9c  e D(Sj'sQ).  This  alternative  characterization  of 

A 

s'."  will  be  subsequently  used  in  the  development  of  the  estimation  procedure. 

^ A n 

Next  let  F(s1,s2, . . . ,sn;  sQ)  = J (Bj(s].;sQ)  + 9cs,.).  By  dropping 

both  the  integrality  and  nonnegativity  restrictions  on  sQ,  we  obtain  the 

following  relaxation  of  Problem  4: 


min  {Qcs  + min  {F(s, ,...,s  is.):  s.  fixed}}, 

s 0 s =0  1 1 n0  0 

j • • • 


If  sQ  is  the  optimal  solution  to  Problem  4,  then 


3F 

8s, 


+ 0c  = 0. 


But 


»r  n 3B.  3T . 

3F  _ y __1 1 

3s  " 3T . 3s 
0 3=1  : 0 


Furthermore,  by  writing  B.(s.;s  ) as  £ Kp(K+s . | X .T . (s  ) ) , 

3 3^  K=1  3 3 3 0 


(8) 


(9) 


we  see  that 


K+s.-l 

3Q.  » -X .T . (s  ) (X.T.(s  ))  3 

— i = y x Ke  33  0 -3-i-0— y 

STj  K=1  1 K=1 


K+s . 

-A.I.(s0,  ».,T.(s B»  1 

3 (K+Sj)! 


= -XjAB(Sj  »sQ). 

Let  B0(sQ)  = y (x-sQ)p(x jXD).  Data  gathered  during  the  conduct  of 


x>s. 


the  study  showed  that  B0(sQ)  can  be  closely  approximated  by  an  exponential 
function.  This  result  should  not  be  entirely  surprising.  The  only  probabilities 
entering  the  backorder  calculation  are  the  so-called  "tail"  probabilities. 

In  the  upper  tail,  the  Poisson  distribution  behaves  almost  like  a geometric 
distribution,  that  is,  each  succeeding  probability  is  approximately  a constant 
proportion  of  the  preceding  value.  This  is  particularly  true  for  large 
values  of  sQ.  Consequently,  when  sQ  is  large  the  expected  number  of 
backorders  existing  at  any  point  in  time  at  the  depot  is  approximately  a 
geometric  function  of  sQ.  Thus  an  exponential  function  provides  a continuous 
approximation  to  the  relationship  between  depot  backorders  and  depot  stock. 

We  therefore  will  approximate  BQ(s0)  ^y  the  exponential  function  of  the  form 


- ~ 1 1 — | 


— ... -A 
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-b  s 

a^e  , where  and  b^  are  positive  real  numbers. 

Then  using  this  approximation  we  rewrite  Tj(sQ)  as 

Tj(s0)=r.Aj.U-r.)(B.  tia0eb°S°) 


and  observe  that 


!!i . . iini . b e-b°s° 

3sq  X To 


Upon  combining  these  observations  we  see  that 


(1-r. ) ~b  s 


^ xj  1 B(sj;V  T1* 


Recall  that  -9c  £ D(s\';s0>.  Consequently  -8c  approximates  the 
marginal  reduction  in  backorders  at  base  j when  the  stock  level  at  that 
base  is  sT.  After  making  this  substitution  and  representing  this  further 

3F  3F 

approximation  of  -r — by  -r — , we  see  that 

dS^  dS  ^ 

o o 

~ = - U 1-r  )X  r 8ca0b0e  - -Oca^e 

0 3=1  J J 

Substituting  this  approximation  into  Equation  9 we  obtain  the  following 
estimate  of  the  optimal  depot  stock  level: 

s = - ~ ln{ — —}  . (10) 

0 bo  aobo 

Recall  the  value  of  sQ  was  derived  based  on  an  exponential  approximation 
of  B^(s^).  As  the  average  number  of  units  in  the  depot  repair  cycle  increases, 
that  is,  as  XD  increases,  the  quality  of  this  exponential  approximation 


improves  in  the  region  in  which  the  optimal  depot  stock  level  should  be  located. 
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Consequently , the  approximation  should  be  most  accurate  in  these  cases.  But 
the  items  for  which  the  search  for  the  optimal  depot  stock  level  is  most  time 
consuming  for  both  the  Sherbrooke  and  Fox-Landi  algorithms  are  the  items 
having  a large  number  of  units  in  depot  repair.  Therefore  the  proposed 
approximation  method  will  be  most  appropriate  for  those  items  requiring  the 
greatest  amount  of  computation. 

The  approach  we  have  described  for  estimating  the  optimal  depot  stock 
level  has  been  coded  and  tested  using  a sample  of  40  items  found  on  the 
Air  Force's  new  F-15  aircraft.  The  test  consisted  of  two  sets  of  runs.  In 
the  first  set,  monthly  flying  was  divided  among  3 bases;  in  the  second  set, 
the  same  monthly  flying  program  was  divided  among  5 bases.  The  total  budget 
distributed  among  the  40  items  ranged  from  $34  million  to  $65  million  in  the 
first  set  of  runs,  and  from  $34  million  to  $88  million  in  the  second  set. 

Table  1 contains  the  data  indicating  both  the  optimal  and  estimated  depot 
stock  level  for  each  item  in  both  runs . 

The  data  show  that  there  usually  is  no  single  optimal  depot  stock  level 
for  an  item.  Rather  the  optimal  value  depends  on  the  amount  of  total  item 
system  stock  available  for  distribution  among  the  depot  and  bases.  The  data 
displayed  in  Table  1 also  indicate  the  estimate  of  optimal  depot  stock  is 
quite  close  to  the  optimal  vaiue  in  all  cases.  Furthermore,  the  increase 
in  expected  system  backorders  using  the  estimated  depot  stock  levels  rather 
than  the  optimal  levels  is  generally  small.  For  most  items  the  increase  is 
substantially  lens  than  .1  backorders. 

The  results  of  the  tests  indicate  that  it  is  possible  to  estimate  closely  the 
optimal  depot  stock  level  using  Equation  10.  Additionally,  incorporating  this 
method  for  estimating  the  optimal  depot  stock  into  the  Fox-Landi  and  Sherbrooke 
algorithms  will  reduce  the  search  required  to  find  the  optimal  depot  stock  level 


i iLrt.  L 
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and  will  therefore  reduce  the  computational  time  needed  to  solve 
Problem  1 using  these  algorithms.  In  the  tests,  the  computation  time  for 
high  demand  items  (XD  >.  20)  was  reduced  by  over  90  per  cent.  Consequently, 
this  approximation  method  is  recommended  for  use  particularly  in  situations 
where  expected  depot  demand  is  large. 


r 
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> . 


Estimated 


Item 

Case  I 
(3  Bases) 

Case  II 
(5  Bases) 

Optimal  Depot 
Stock  Levels 

1 

4-7 

5-9 

6 

2 

1,2 

1-3 

1 

3 

6 

6,7 

6 

4 

0-2 

2,3 

1 

5 

10,11 

8-12 

10 

6 

18-21 

18-21,25 

19 

7 

1,2 

1,2 

1 

8 

2 

3,4 

2 

9 

5,6 

6,7 

6 

10 

1 

1,2 

1 

11 

4,5 

4-6 

5 

12 

1 

1 

0 

13 

0-2 

0,1 

0 

14 

1-3 

1-3 

2 

15 

2-4 

3,4 

3 

16 

8,9 

8,9 

8 

17 

1,2 

1,2 

1 

18 

3,4 

3-5 

3 

19 

12-14 

13-14 

12 

20 

9-12 

10-13 

10 

21 

21-27 

22-28 

23 

22 

4,5 

4-6 

5 

23 

1 

1-3 

1 

24 

1,2 

2,3 

2 

25 

5-7 

6,7 

6 

26 

16 

16 

16 

27 

3 

3,4 

3 

28 

40-42 

41-43 

40 

29 

8-10 

9,10 

9 

30 

1 

2 

9 

31 

1,2 

1,2 

1 

32 

8,9 

8,9 

8 

33 

4,5 

5,6 

5 

34 

9-11 

9,10 

10 

35 

6,7 

7,8 

7 

36 

1-3 

2 

2 

37 

1,2 

1,2 

1 

38 

7,8 

7-9 

7 

39 

2,3 

3,4 

3 

40 

41-43 

42-44 

41 

TABLE  1 


Comparison  of  Optimal  and  Estimated 
Depot  Stock  Levels 
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n 1968  Sherbrooke  formulated  the  well  known  METRIC  model  for  a two-echelon 
inventory  system  consisting  of  a set  of  bases  and  a supporting  depot.  The 
items  stocked  in  the  system  are  called  recoverable  items,  that  is,  they  are 
subject  to  repair  when  they  fail.  Sherbrooke  and  others  have  proposed  a variety 
of  algorithms  for  determining  optimal  base  and  depot  stock  levels.  A substantial 
portion  of  the  computational  requirement  associated  with  each  of  the  algorithms  — j 
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is  related  to  the  search  for  the  optimal  depot  stock  level.  The  purpose 
of  this  paper  is  to  describe  an  easily  implementable  method  for  estimating 
the  optimal  depot  stock  level.  The  computational  experience  reported  in  the 
paper  indicates  the  proposed  method-provides  an  excellent  estimate  of  the 
optimal  depot  stock  level  particularly  for  high  demand  items,  furthermore, 
the  proposed  method .significantly  reduces  the  computational  requirements  for  any 
known  algorithm  for  solving  Sherbrooke's  problem. 
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